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' The QCD vacuum can be studied using external fields. We report here results respectively obtained 

^— «, ' probing the lattice QCD vacuum by means of an abelian monopole field and of an abelian chromo- 

' magnetic field. 
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1 Introduction 



Understanding color confinement is one of the long standing problems of the Standard 
Model of Particle Physics. Notwithstanding many theoretical efforts a totally convincing 
explanation of the confinement phenomenon is still lacking (we refer the reader to recent 

reviews on confinement ^ ' ^ ' ^ ' ^ j) and a full understanding of the QCD vacuum dynamics 
^ . is not yet at our disposal. Nevertheless lattice computations can shed light on many 

aspects of the QCD vacuum and possibly help in identifying the mechanism of the color 
^ , confinement: that has surely been a leading interest in Adriano's research activity. In 

■ this paper we review some results we have collected in the last few years in studying the 

Q I QCD vacuum by means of external fields. The plan of the paper is the following. In 

^ ■ section|2lwe briefly recall our gauge invariant method to put external fields on the lattice. 

O I In section |3] we consider the QCD vacuum under the influence of an Abelian monopole 

external field and we study the phase transition, even with two dynamical flavors, using 
the free energy in presence of the monopole field. In section |4] we probe the vacuum 
structure by means of a constant abelian chromomagnetic field, both at zero and finite 
temperature. Finally in section|5lwe present our conclusions. 



2 External Fields on the Lattice 

2. 1 Pure gauge theories at zero temperature. The lattice effective action 

In previous papers l^^^ we introduced a lattice effective action r[74^'"] for an external 
background field A^^\ which is gauge invariant against static gauge transformations of 
the background field 

r[i'-t] = -lln/^^| (1) 



m j 

where Lj is the lattice size in time direction and A^^\x) is the continuum gauge potential 
of the external static background field. is the lattice partition functional 

ZlA^""^] = I VU e-^^ , (2) 

with Sw the standard pure gauge Wilson action. 



The functional integration is performed over the lattice links, but constraining the 
spatial links belonging to a given time slice (say xt = 0) to be 

Uk{x, xt = 0) = Ul^\x) , {k = l, 2, 3) , (3) 

[/^"'(x) being the lattice version of the external continuum gauge potential A^'^^x) = 

A^^^{x)Xa/2. Note that the temporal links are not constrained. 

In the case of a static background field which does not vanish at infinity we must also 
impose that, for each time slice xt 7^ 0, spatial links exiting from sites belonging to the 
spatial boundaries are fixed according to eq. Q. In the continuum this last condition 
amounts to the requirement that fluctuations over the background field vanish at infinity. 

The partition function defined in eq. ^ is also known as lattice Schrodinger func- 
tional and in the continuum corresponds to the Feynman kernel^. Note that, at variance 

with the usual formulation of the lattice Schrodinger functional^ where a lattice cylindri- 
cal geometry is adopted, our lattice has an hypertoroidal geometry so that Sw in eq. Q 
is allowed to be the standard Wilson action. 

The lattice effective action r[74'^'"] corresponds to the vacuum energy, ii^o[^'^'"]> in 
presence of the background field with respect to the vacuum energy, £^0 [0] , with A^^^ = 0. 



2.2 Pure gauge theories at finite temperature. The thermal partition fiinctional 

If we now consider the gauge theory at finite temperature T = 1/ {aLt) in presence of an 
external background field, the relevant quantity is the free energy functional defined as 



1 J 



T[A<^'''] = --\n{ r> • (4) 
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Zt [A^^^] is the thermal partition functionally in presence of the background field A^^^, 
and is defined as 



A^""' = VU e"^'^ . (5) 

J JUk{xM)=Uk{xfi)=W^\x) 

In eq. as in eq. the spatial links belonging to the time slice = are constrained 
to the value of the external background field, the temporal links are not constrained. The 
free energy functional eq. @ corresponds to the free energy, in presence of the 

external background field evaluated with respect to the free energy, -F[0], with A'^^^ = 0. 
When the physical temperature is sent to zero the free energy functional reduces to the 
vacuum energy functional eq. (0. 

2.3 Including dynamical fermions 

When including dynamical fermions, the thermal partition functional in presence of a 
static external background gauge field, Eq. (|5Jl, becomes: 



Zt 



A^ 



t/fc{Lt,a-)=C/fc(0,a;)=C/f (a;) 

VUe-^^ det M , (6) 

Uk{Lt,x)=Uk{0,x)=Ufix) 



where Sw is the Wilson action, Sp is the fermionic action and M is the fermionic ma- 
trix. Notice that the fermionic fields are not constrained and the integration constraint 
is only relative to the gauge fields: this leads, as in the usual QCD partition function, 
to the appearance of the gauge invariant fermionic determinant after integration on the 
fermionic fields. As usual we impose on fermionic fields periodic boundary conditions 
in the spatial directions and antiperiodic boundary conditions in the temporal direction. 



3 Abelian Monopoles 

A mechanism for color confinement based on dual superconductivity of the QCD vacuum 

by abelian monopole condensation has been proposed a long time ago 10 11 ^ disorder 
parameter which is related to abelian monopole condensation in the dual superconduc- 
tivity picture of confinement has been developed by the Pisa group and consists in the 
vacuum expectation value of a magnetically charged operator, (fi) . 

Our pro posal is to detect abelian monopole condensation by looking at the free en- 

ergyll^ll^in presence of an abelian monopole background field. 

Since the free energy is related to the vacuum dynamics and not, like the trace of the 
Polyakov loop, to a symmetry of the gauge action which is washed out by the presence 
of dynamical fermions, we feel that it could be used to detect the finite temperature phase 
transition also in full QCD. 



3. 1 Abelian monopole condensation in pure gauge theories 

For SU(3) gauge theory the maximal abelian group is U(l)xU(l), therefore we may 
introduce two independent types of abelian monopoles using respectively the Gell-Mann 
matrices A3 and Ag or their linear combinations. 

In the following we shall consider the abelian monopole field related to the A3 diago- 
nal generator. In the continuum the abelian monopole field is given by 



2 X ( X — X ■ n) 



where n is the direction of the Dirac string and, according to the Dirac quantization 
condition, nmon is an integer. The spatial lattice links are constrained (see Eq. ^) to the 
values corresponding to the continuum abelian monopole field Eq. Q. 

The monopole background field is introduced by constraining the spatial links exiting 
from the sites at the boundary of the time slice = (see Eq. For what concern 
spatial links exiting from sites at the boundary of other time slices {xt 7^ 0), since the 
monopole field vanishes at infinity, we consider two possibilities. In the first one we 
constrain these links according to Eq. Q (in the following we refer to this possibility as 
"fixed boundary conditions"). In the second one we do not impose the constraint Eq. Q 
on the above mentioned links (this possibility will be referred as "periodic boundary 
conditions"). 

We simulate pure SU(3) lattice gauge theory with Wilson action. The lattice geometry 
is hypertoroidal. In Fig. ^ we report our numerical results for the derivative of the free 
energy F' . The data display a sharp peak that increases by increasing the lattice spatial 
volume. By inspecting Figure[2it is evident that F'{f5) = in a finite range of /?, starting 
from /? = and below the critical coupling, signaled by a peak in F'{f3). Therefore 
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Figure 1: F' vs. /3 (Wmon = 10) for Ls = 16, 24, 32 and Lt = 4. 



F{(3) = in a finite range below the critical coupling, above which the gauge system gets 
deconfined. The vanishing of the free energy implies abelian monopole condensation. 
On the other hand, F'{(3) becomes different from zero and increases with the lattice 
spatial volume near the critical coupling, as expected in presence of a phase transition. 
Moreover in the weak coupling regime F'{(3) stays constant and almost independent 
from the spatial lattice volume. This corresponds to the classical monopole energy which 
depends linearly on /?. It is clear that in the deconfined phase it costs a finite amount of 
energy to create a monopole, and for, there is not abelian monopole condensation. 

A finite size scaling analysis can be performed and it can be shown that the scaling 
parameters are compatible with a first order phase transition. As a first step we determine 
the value of the critical coupling Pc{Lf^), where Lf^ is the effective spatial size (Lf^ = 

— 2 in the case of fixed boundary conditions). We fitted our lattice data with the scaling 
law 

r^'fp reffN _ "■ij^f)'^ „ 



The output of the fits are reported in Table Q The parameter di has been fixed at the 




value obtained with the fit Eq. We see clearly that a/u agrees with 7 within statistical 
errors, as expected if we want a sensible result in the thermodynamical limit Lf^ 00 
(see Eq. 

In Fig. 121 we plot F'{/3,Lf^) rescaled with the factor l/{Lf^y versus the scaling 
variable x, respectively for spatial "fixed boundary conditions" and for spatial "periodic 
boundary conditions". The full line is 

F'{P,Lf) ^ ai 

We find that the scaling relation holds quite well for a very large range of x. The quality 
of the scaling can be inferred looking at Fig. |2l 

It is interesting to comment on the behavior of ex p(— F(/3)/ T), which is the analo- 
gous of disorder parameter developed by the Pisa groupP^l^^l in the thermodynamical 
limit implied by Eq. ©. Indeed we have that 

exp (-^) = exp (^-^ 1^ F'iP') dp'^ , /3o < /? < /3c , (10) 



spatial "fixed boundary conditions" 




I3c v di 


a 


12.199 


1.247 5.3251 0.335 0.6 


0.351 


±3.9004 


±0.089 ±0.0110 ±0.026 constant 


±0.035 


spatial "periodic boundary conditions" 
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-1 j3c y di 


a 


13.461 


1.510 5.3222 0.340 0.6 


0.347 


±1.337 


±0.555 ±0.0013 ±0.020 constant 


±0.009 



Table 1 : The values of the parameters obtained by fitting Eq. ^8} to the data for the derivative of the monopole 
free energy on lattices with spatial volumes 16'', 24^, and 32'' and spatial "fixed" or "periodic" boundary 

conditions respectively. 



while we already know that exp(— F(/3)/T) 
size (see Fig. Q). From Eq. (ISj we get 



Ifor (3 < irrespective of the lattice 



exp 



FiP) 



exp 



1 

T 



Ol 



1 



a 



{\/3o-Pc 



ll-a 



l/?-/3c|'-") 



l3o<l3<Pc 



So that F{f3) decreases when 
On the other hand, for a = 1 
power of (/5c — /3). Thus we 
exp{-F{(3)/T) at (3c and the 
ever, it must be stressed that 
small so that exp{-F{/3)/T) 



(11) 

/3— >/3cif0<a<l tending to a finite value at (3 = f3c- 
it is easy to see that exp(— F(/3) /T) decreases to zero as a 
conclude that for a < 1 we have a discontinuous jump of 
strength of the discontinuity weakens when a I. How- 
the discontinuous jump of exp(— F(/3)/T) is exceedingly 
decreases almost continuously toward zero when /3 — > f3c- 



3.2 The deconfinement transition in QCD with two dynamical flavors 

We compute the derivative of the monopole background field free energy with respect 
to the gauge coupling where now the expectation value is evaluated with the full QCD 
action. Our main goal is to try to use our data for F'{(3) on different spatial volumes and 
different bare quark masses to infer the critical behavior of two flavors full QCD near the 

deconfining transition (for details on numerical simulations see Ref.^l). 

We perform a finite size scaling analysis using the following scaling law (suggested 
by Eq. (Hi 

\{Lfy/-^{P-Pc{m,))-d,\- ^''^ 

where the critical coupling I3c{mq) depends on the quark mass niq. The dependence 
of the critical coupling j3c{mq) on the quark mass is determined by the chiral critical 

point 19 20 j-j^g thermodynamical limit, by known universality arguments the criti- 
cal couplings will scale like 

f3c{mq) = (3c{mq = 0) ± cm]"^^ , (13) 




Figure 3: F'{(3, if, m,) rescaled by the factor (Lf )^. The values of if and niq — am are displayed in 

the legend. 



where is a combination of critical exponents which for the case of 2-flavors QCD 
are expected to be those of the three-dimensional 0(4) symmetric spin models l//3(5 ~ 
0.5415. Inserting Eq. (fT3l into Eq. (fT2l we are lead to the following scaling law 

^ ' " |(Le/)i/^(/?-/3,(0)-cm^)-dir ' ^ ^ 

where again j = a/u assures a sensible thermodynamical limit. Note that, to take care 
of finite volume effects, the exponent rj is expected to be = i^ c/i^ with Uc = i^'/PS, 
where u', (i, and 6 are the chiral critical exponents. Indeed, Eqs. (fl4t assures that in the 
scaling region 

""^^^Lf'^'^ = ^^^^f)""'^^ - PM), {Lf)'''^m,) . (15) 

In our case the relevant chiral critical exponents are those of the three-dimensional 0(4) 

symmetric spin models where v' = 0.7423 , Uc = 0.4019. In Table|2lwe report the 
results obtained by fitting Eq. (I14t to all our lattice data. 

From Table |2l we can see that a/u = 2.35 it 0.34 consistent with 7 = 2.00 it 0.47. 
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79.4 


2.00 


4.9933 


0.54 1.10 0.31 


0.6 


0.728 


±76.6 


±0.47 


±0.0138 


±0.11 ±0.19 ±0.03 


constant 


±0.078 



Table 2: The values of the parameters obtained by fitting Eq. <14> to the data for the derivative of the 
monopole free energy in two-flavors full QCD on lattices with spatial volumes 16'', 24"', and 32^ and Lt — 4. 



Concerning the parameter rj we find that it is poorly determined by our data. If we con- 
strain 7] in our fit we get rj = 1.10 ± 0.19 which, together with u = 0.31 ± 0.03 leads 
to Uc = 0.34 ± 0.07 consistent with the value reported above. However if we release the 
constraint on rj our data can also be fitted with smaller values for rj without altering sig- 
nificantly the other parameters. Moreover by confronting the exponent a in Table|2lwith 
the corresponding value for the SU(3) pure gauge in Tabled we conclude that our data 
for full QCD with two dynamical flavors are compatible with a first order phase transition 
= 0.31 ± 0.03) but this is weaker than in the quenched case. Our results are in agree- 
ment with the indications for a first order phase transition in ful l QCD with 2 dy namical 

flavours (in the same range of quark masses) obtained elsewhen 
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4 Abelian Chromomagnetic Field 

We feel that it could be useful to study vacuum dynamics by using different perspec- 
tives in order to shed light on the basic mechanism of color confinement. Indeed, as 

recently observed^ in connection with dual superconductivity picture, even if magnetic 
monopoles do condense in the confinement mode, the actual mechanism of confinement 
could depend on additional dynamical forces. 

We found ^ that a weak constant abelian chromomagnetic field at zero temperature is 

completely screened in the continuum limit, while at finite temperature!^ our numerical 
results indicate that the applied field is restored by increasing the temperature. These 
results strongly suggested that the confinement dynamics could be intimately related to 
abelian chromomagnetic gauge configurations. 

We report here numerical results showing that for non abelian gauge theories the de- 
confinement temperature depends on the strength of an applied external constant abelian 

chromomagnetic fieldESI This is at variance of abelian magnetic monopoles where the 

abelian monopole background fields do not modify the deconfinement temperature^^. 

4. 1 The color vacuum Meissner effect 

Let us now define a static constant abelian chromomagnetic field on the lattice. We 
consider the SU(3) case in (3+1) dimensions, a thorough discussion of the SU(2) and 

U(l) case can be found elsewhere . 

In the continuum the gauge potential giving rise to a static constant abelian chromo- 
magnetic field directed along spatial direction 3 and direction a in the color space is given 
by 

A'^\x)=A'^\x)5a,a, At\x) = 5k,2XlH . (16) 



The constrained lattice links (see Eq. ^) are obtained accordingly to the continuum 
gauge potential Eq. (fT6t . Since our lattice has the topology of a torus, the magnetic field 
turns out to be quantized 
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next integer. 



(17) 



In the case of a constant background field the relevant quantity is the density /[A*^"'] 
of free energy. We evaluate by numerical simulations the derivative with respect to the 

coupling f3 of the free energy density / [A^^^] at fixed external field strength gH. 

We consider here SU(3) pure gauge theory. As is well known, the pure SU(3) gauge 
system undergoes a deconfinement phase transition at a given critical temperature. The 

critical coupling /3c can be evaluated by looking at the peak of /'[A*^"'], the derivative of 
the free energy density with respect to /3. 

Once the critical coupling P*{Lt) has been determined, the deconfinement tempera- 
ture can be given in units of the string tension Tc/ \J cr{f3c) = 1/Lt \/a{l3c)- Moreover, 
using eq. (fTTb . the field strength is ^ gH j cr{(3c) = \/47rnext/ixO'(/9c)- Our data for 
Tc/ versus \/gH/ on a 64'^ x 8 lattice are displayed in Fig. 0] It is worth to note 
that lattice data can be reproduced by the linear fit 



with Tc{0)/./a = 0.643(15) and a = -0.245(9). The critical field can now be ex- 
pressed in units of the string tension and assuming ^/a = 420 MeV, the critical field is 
^/gH^ = (1.104 ± 0.063)GeV corresponding to gHc = 6.26(2) x 10^^ Gauss. 

Therefore we may conclude that a critical field exists such that Tc = for gH > gHc- 
This kind of behavior could be interpreted as the colored counterpart of the Meissner 
effect in ordinary superconductors, when strong enough magnetic fields destroy the su- 

perconductive BCS vacuum . Then we shall refer to this remarkable result as the 
reversible color Meissner effect. It is worthwhile to stress that this effect is not related 
to the color superconductivity in cold dense quark matter. Indeed, we believe that our 
reversible color Meissner effect is deeply rooted in the non-perturbative color confining 
nature of the vacuum and could be a window open towards unraveling the true nature of 
the confining vacuum. 

5 Conclusions 

We investigated the QCD vacuum using external fields. 

We first considered the nature of deconfining phase transition in SU(3) pure gauge the- 
ory and full QCD with two flavors of staggered fermions. To locate the phase transition 
we used the derivative of the monopole free energy with respect to the gauge coupling. 
The monopole free energy is defined by means of a gauge invariant thermal partition 
functional in presence of the abelian monopole backgroud field. In the pure gauge case 
our finite size scaling analysis indicate a weak first order phase transition. In the case 
of 2 flavors full QCD, we performed simulations by varying spatial lattice sizes and 
quark masses. We find that deconfinement transition in full QCD with 2 degenerate 




(18) 



Figure 4: SU(3) in (3+1) dimensions. The critical temperature Tc estimated on a 64^ x 8 lattice in units of 
the string tension versus the square root of the field strength \/ gH in units of the string tension. Solid line 
is the linear fit eq. <18> . In correspondence of zero vertical axis: open circle is Tc/ at zero external field; 
full circle is the determination of Tc/^/o obtained in the literature. 



dynamical flavors is consistent with a weak first order phase transition, contrary to the 
expectation of a crossover for not too large quark masses, but in agreement with recent 

indications^^ 

We then focused our attention to (3+1) dimensional SU(3) pure gauge theory in pres- 
ence of an uniform chromomagnetic field. We reported numerical evidences in favor of 
what we named "color vacuum Meissner effect" that implies the existence of a critical 
field gHc such that for gH > gH^ the gauge system is the deconfined phase. Such an ef- 
fect seems to be generic for non abelian gauge theories^^ (see Fig.|5Jl. We speculate that 
the peculiar dependence of the deconfinement temperature on the strength of the abelian 
chromomagnetic field gH could be naturally explained if the vacuum behaved as an 
ordinary relativistic color superconductor, namely a condensate of color charged scalar 
fields whose mass is proportional to the inverse of the magnetic length. However, the 
chromomagnetic condensate cannot be uniform due to gauge invariance of the vacuum, 
which disorders the gauge system in such a way that there are not long range correlations. 
Consequently if the vacuum behaved as a non uniform chromomagnetic condensate, our 
reversible color Meissner effect could be easily explained, for strong enough chromo- 
magnetic fields would force long range color correlations such that the gauge system gets 
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Figure 5: Phase diagram of four dimensional SU(2) and SU(3) gauge theories. 

deconfined. One might thus imagine the confining vacuum in non abelian gauge systems 
as a disordered chromomagnetic condensate which confines color charges due both to 
the presence of a mass gap and the absence of long range color correlations, as argued by 

an 

R.P. Feynman in (2+1) dimensions^^. 
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